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Abstract. We present a new implementation of the magnetic Doppler imag- 
ing technique, which aims at self-consistent temperature and magnetic mapping 
of the surface structures in cool active stars. Our magnetic imaging procedure 
is unique in its capability to model individual spectral features in all four Stokes 
parameters. We discuss performance and intrinsic limitations of the new mag- 
netic Doppler imaging method. A special emphasis is given to the simultaneous 
modelling of the magnetically sensitive lines in the optical and infrared regions 
and to combining information from both atomic and molecular spectral features. 
These two techniques may, for the first time, give us a tool to study magnetic 
fields in the starspot interiors. 



1. Introduction 

Rotational modulation of the intensity and polarization spectra of active stars 
offers a unique possibility to reconstruct 2-D maps of the photospheric magnetic 
fields and temperature spots, thus providing a key constraint for the theoret- 
ical modelling of the stellar activity an d dynamo. Until now Zeeman Doppler 
imaging of the ac tive late-type stars (iBrown et al.l Il99ll ; iHussain et al. 2000; 



Donati et al.l 120031 ) was based upon a non-simultaneous and greatly simplified 



interpretation of the average (least-squares deconvolved) Stokes / and Stokes 
V spectra. In contrast, detailed and physically c onsistent methods of the stel- 
lar magnetic field mapp i ng ha ve been developed ( Piskunov fc Kochukhovl [20021 ; 



Kochukhov fc Piskunovl 120021 ) and successfully applied to the high-resolution 



circular and linear polarizati on observations of the early-type magnetic stars 
flKochukhov et al]l200ll2004h . 



Improvements in the observational techniques and instrumentation have re- 
cently made it possible to de tect polarization signat ures inside individual atomic 
( Petitl [20061 ) and molecular ( Berdyugina et a l. 2006) lines of the late- type stars. 



With these observational data in mind, we have developed a new Magnetic 
Doppler imaging (MDI) code Inversl3. This novel tool aims at self-consistent 
reconstruction of the temperature and magnetic structures in cool active stars. 
A special property of our magnetic DI procedure is its capability to realisti- 
cally model individual atomic and molecular lines in all four Stokes parameters 
simultaneously. Here we describe physical foundations and numerical implemen- 
tation of our MDI technique. We also present numerical experiments designed 
to evaluate performance of Inversl3. 
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2. Magnetic Doppler imaging with Inversl3 

The new MDI code reconstructs iteratively two-dimensional distributions of tem- 
perature and three components of the magnetic field vector using high-resolution 
spectropolarimetric time series recorded in the Stokes / and V or all four Stokes 
parameters. Inversions employ the Levenberg-Marquardt minimization proce- 



dure and are constrained with the Tikhonov regularization (jPiskunov fc Kochukhov 



2002). 

Our modelling of the magnetic field and cool spots is self-consistent, i.e. 
the local intensity and polarization profiles are computed taking into account 
the Zeeman effect and the presence of temperature inhomogeneities at the same 
time. This allows us to deduce correct field strength inside cool spots. In 
contrast, all previous ZDI methods, that use only the Stokes V for magnetic 
mapping, significantly underestimate magnetic flux in cool spots or fail to detect 
it altogether. 

The magnetic spectrum synthesis is carried out using realistic model at- 
mosphere grids of cool stars (ATLAS9, Phoenix, MARCS). The atomic and 
molecu lar concentrations are calculated with the advanced equation of state 
solver (IValenti et al.lll998l ). The local Stokes IQUV parameters are computed 



by solvin g the polarized radiative t ransfer equation with the qua dratic DELO 
method ( Socas-Navarro et a l. 2000; [pTskunov fe Kochukhov! [2002). 



3. Numerical simulations 

Performance of the new MDI code is verified with the numerical experiments. In 
these tests Inversl3 infers temperature and magnetic maps from the simulated 
observational data in two or four Stokes parameters. For all tests presented here 
we used the stellar parameters T e ^ = 5000 K, logg = 4.0, vsini = 40 kms -1 , 
i = 50°. We assumed the temperature contrast of 1500 and 750 K for the umbral 
and penumbral starspot regions, respectively. The magnetic field strength was 
set to 2 kG in the center of a spot, reducing to 1 kG in the penumbra. Simulations 
were performed for the three optical magnetically sensitive Fe I lines at A 5497.5, 
5501.5, 5506.8 A, two IR lines Fe I 15648.5, 15652.9 A, and the TiO bandhead 
at 7055 A. 

In the first simulation described here we studied the capability of Inversl3 
to infer temperature contrast of cool starspots. Fig.[T]shows reconstruction of the 
temperature distribution consisting of four spots located at different latitudes. 
For these large-scale features the spot-to-photosphere temperature contrast is 
underestimated by no more than 200-300 K even when we use only the Fe I lines 
around 5000 A. The temperature reconstruction is improved (errors reduced to 
50-100 K) if the TiO band or the IR lines are included in the inversion. 

In the second series of tests we have performed magnetic inversions for large, 
purely magnetic spots containing radial, meridional and azimuthal field. Inver- 
sion for the case of radial magnetic topology is illustrated in Fig. El It is evident 
that the radial field is successfully recovered for the high- and intermediate- 
latitude spots. At the same time, a strong radial to meridional field crosstalk 
occurs for the low-latitude spots in the Stokes IV imaging. This problem is 
alleviated when linear polarization data are included in the magnetic inversion. 



<p=0.00 



Magnetic Doppler Imaging of Active Stars 
<p=020 (p=OAO <p=0.60 (p=O.SO 




Figure 1. Reconstruction of the temperature inhomogeneities with 
Inversl3. Upper row: the input temperature distribution. Lower row: re- 
constructed map. 
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Figure 2. Reconstruction of the magnetic field topology consisting of four 
spots with radial field. Upper row: the true distribution of magnetic field. 
Middle row: reconstruction using the Stokes IV spectra. Lower row: recon- 
struction using all four Stokes parameters. 



Finally, we have tested possibility to probe magnetic fields inside cool 
starspots. In Fig. [3] we illustrate the adopted structure of the temperature 
and magnetic inhomogeneities and show results of different magnetic inversions. 
We find that if the magnetic field is studied using the Stokes V alone, without 
accounting for the temperature inhomogeneities (as it is done in all previous 
implementations of ZDI), the resulting map misses entirely the magnetic flux 
in the regions of reduced temperature. On the other hand, our standard MDI 
reconstruction based on the self-consistent interpretation of the Stokes / and V 
spectra recovers correctly the geometry of the field, although the field strength 
is underestimated by about 50% in the starspot center. The agreement between 
the true and reconstructed maps can be improved significantly if the IR Fe I 
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Figure 3. Reconstruction of the magnetic field inside and outside cool 
starspots. The 2-D images show distribution of the radial magnetic com- 
ponent, whereas the contour lines show location of cool spots. The spot at 
the longitude 280° combines reduced temperature and a strong magnetic field. 
The other two spots are magnetic or temperature only. Different rectangular 
plots show the true maps of surface inhomogeneities (upper row) and the re- 
constructions (lower row) using Stokes V alone (left column), self-consistent 
Stokes IV mapping (middle column) with the optical lines, and a simultane- 
ous IV mapping with the optical and IR lines (right column). 



lines are incorporated in the inversion. This becomes possible due to a much 
smaller spot-to-photosphere flux contrast in the IR and the resulting larger rela- 
tive contribution of the spot polarization signatures to the disk-integrated Stokes 
V spectra. The polarization diagnostic of the TiO molecular bands allows us to 
achieve similar improvement in the magnetic mapping of the starspot regions, 
although the MDI with molecular lines suffers from reduced spatial resolution 
due to a large width of molecular features. However, the main advantage of 
the simultaneous analysis of molecular and atomic lines is the possibility to use 
the same observational material for the inversion, whereas costly multisite cam- 
paigns are needed to secure simultaneous visual and IR high-resolution spectra. 
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